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Boron clusters have been widely studied theoretically for their geometrical properties and elec- 
tronic structure using a variety of methodologies. An important cluster of boron is the B12 cluster 
whose two main isomers have distinct geometries, namely, icosahedral (Ih) and quasi planar (Cz v ). 
In this paper we investigate the linear optical absorption spectrum of these two B12 structures with 
the aim of examining the role of geometry on the optical properties of clusters. The optical absorp- 
tion calculations are performed using both the semi-empirical and the ab initio approaches. The 
semi-empirical approach uses a wave function methodology employing the INDO model Hamil- 
tonian, coupled with large-scale configuration interaction (CI) calculations, to account for the 
electron-correlation effects. The ab initio calculations are performed within a time-dependent- 
density-functional-theory (TDDFT) methodology. The results for the two approaches are in very 
good qualitative agreement with each other. Quantitatively speaking, results agree with each other 
in the lower energy region, while in the higher energy region, features predicted by the TDDFT 
approach are red-shifted as compared to the INDO-CI results. Both the approaches predict that 
the optical absorption begins at much lower energies in the icosahedral cluster as compared to the 
planar one, a fact which can be utilized in experiments to distinguish between the two geometries. 
At higher energies, both the isomers exhibit plasmon-like excitations. 

PACS numbers: 36.40.Vz, 31.10.+z,31.15.bu,31.15.vq 



I. INTRODUCTION 

In the recent years, the study of the boron clusters 
has attracted considerable attention due to their po- 
tential applications as hydrogen storage devices hard 
semiconducting solids and various other properties 
Extensive theoretical studies of one-, two-, and three- 
dimensional boron clusters of diff erent sizes have been 
carried out by various researchers! 5 | 6 | 7 | 8 | 9 | 1Q | 11 1 12 1 13 1 14 1 15 1 
Of these, icosahedral B i2 cluster^ has generated a lot 
of interest in recent years possibly because of its high 
symmetry and its occurrence as the fundamental struc- 
tural unit in crystalline boron, and various boron-rich 
solids P^E^kn an early work Bambakidis and Wagner^ 
computed the structural and cohesive properties of icosa- 
hedral B12 using the SCF-Xa-SW approach, and esti- 
mated its HOMO-LUMO gap. Kawai and WarcP em- 
ployed Car-Parinello ab initio molecular dynamics ap- 
proach to analyze the structural instabilities of cage-like 
B12. Kato and Yamashita 7 performed ab initio Hartree- 
Fock calculations to optimize the geometries of vari- 
ous neutral and cat ionic boron clusters, and obtained 
a triplet ground state with a trigonal bipyramid struc- 
ture for the B12 cluster. BoustaniJ^ using both ab ini- 
tio density-functional theory (DFT), 8 and wave function 
based quantum-chemical 9 approaches, presented detailed 
studies of various boron clusters, including several struc- 
tural isomers of B12. Using a semiempirical approach, 
Fujimori and KimursP^ explored the nature of bonding 
in icosahedral clusters of group III atoms, including cage- 
like Bi 2 . Hayami, 11 using an a b initio DFT based ap- 
proach studied the encapsulating properties of the icosa- 
hedral Bi 2 . Zhai et al., 12 in a recent joint theoretical and 
experimental study investigated some small boron clus- 
ters, including B12, and argued that these clusters prefer 



planar aromatic structures. In a first-principles study, 
He et al. explored the ionicities of boron-boron bond 
in icosahedral B12 caused by symmetry breaking. Atis 
et al. based upon an ab initio DFT methodology, pre- 
sented a theoretical study of several neutral boron clus- 
ters, along with an analysis of relative stability of various 
isomers, including B12 isomers. Prasad and Jemmis^in 
a first principles study examined the structure of several 
large boron clusters and concluded that the ones con- 
structed from the B i2 icosahedron as the basic unit are 
more stable as compared to the fullerene-like structures. 

In spite of several studies of structural properties of 
boron clusters, very few studies of their dielectri c re- 
sponse properties exist. For example, Reis et al. j 1 8 1 19 1 
computed the static linear and nonlinear optical suscep- 
tibilities of rhombic B4 clusters using a first-principles 
methodology. In an earlier work, by means of ab ini- 
tio correlated calculations, we had computed the static 
dipole polarizabilities of ladder-like boron clusters.^ Al- 
though optical absorption spectra of icosahedral clus- 
ters of atoms such as Al and Pb have recently been 
computed p2 to the best of our knowledge, no such cal- 
culations for the linear or non-linear optical response 
of B12 clusters (quasi-planar or icosahedral) have been 
performed. Given the fact that B i2 clusters were re- 
cently discovered experimentally,^ it is of considerable 
interest to compute their optical properties. Theoretical 
calculations of optical absorption spectra, coupled with 
the experimental measurements, can be used to identify 
clusters of various shapes and sizes, and to distinguish 
between various isomers. With this aim in mind, we 
present a theoretical study of the frequency-dependent 
linear optical response of B i2 icosahedral cluster, as well 
as its quasi-planar isomer which was recently discovered 
experimentally by Zhai et al\^ Such comparative stud- 
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ies are important because they also help us in under- 
standing the role of geometry in the optical response 
of clusters. Our calculations have been performed us- 
ing both a semiempirical INDO 22 model based approach 
within the configuration-interaction (CI) framework, as 
well as an ab initio time-dependent dens ity functional 
theory (TDDFT) based methodology) 2 ^ 2 -^ Linear optical 
absorption spectra computed using the two approaches 
are in full qualitative agreement with each other. Re- 
garding quantitative aspects, agreement in the lower en- 
ergy region is very good, while the TDDFT results are 
redshifted as compared to the INDO-CI results in the 
high energy region. The predictions of our calculations 
can be tested in future optical absorption experiments on 
these clusters. 

Remainder of the paper is organized as follows. In 
the next section we present the theoretical aspects of our 
approach, to be followed by the presentation and discus- 
sion of our results in section [TTTj Finally, in section [IV] we 
present the conclusions 



II. THEORETICAL BACKGROUND 

For the present study, we adopted a wave- 
function based electron-correlated methodology employ- 
ing the semi-empirical intermediate-neglect of differen- 
tial overlap (INDO) model developed by Pople and 
coworkers! 22 * 25 ^ The INDO model, like its predecessor 
complete neglect of differential overlap (CNDO) modelp^ 
employs an effective valence-electron Hamiltonian which 
uses Slater-type-orbital (STO) as basis functions, and 
several of its one- and two-electron integrals are approx- 
imated using a semi-empirical parameterization scheme. 
The main difference between the INDO and the CNDO 
models is that, in the INDO model one-center Coulomb 
repulsion integrals include more nonzero terms as com- 
pared to the CNDO model, leading generally to a better 
description of the excited state energies. For a detailed 
discussion of the theory behind these approaches, and es- 
sential differences therein, we refer the reader to the book 
by Pople and BeveridgePSl 

We adopted a semi-empirical approach, as against a 
fully ab intio one, because, with twelve atoms, an ab ini- 
tio correlated calculation with even modest-sized basis 
sets becomes intractable. On the other hand, the INDO 
method, with its smaller basis set (four basis functions 
per atom), allows one to treat electron-correlation effects 
at a much higher level than what is feasible using an 
ab initio approach. Our calculations are initiated at the 
Hartree-Fock (HF) level, within the INDO model, us- 
ing a computer program developed recently by us The 
INDO-HF molecular orbitals (MOs), are used to trans- 
form the Hamiltonian from the original atomic-orbital 
(AO) to the MO representation, which is subsequently 
used in the post-HF correlated calculations. 

There are several variants of the INDO approach in 
vogue which differ from each other in terms of the semi- 



empirical parameters used. In this work we have used 
the original INDO parameterization proposed by Pople 
and coworkers. 22 For the calculations of spectroscopic 
properties such as the e xcitation energies, CNDO/£p3 
and the INDO/S- 28 | 29 | 3Q | approaches have been used so 
frequently that it is virtually impossible to cite all of 
them. INDO /S met hod was parameterized by Zerner 
and coworker d 28 * 29 * 30 * to reproduce spectra of small aro- 
matic molecules with the use of low-order CI approaches, 
such as the singles-CI (SCI) method. However, when 
Adachi and Nakamura benchmarked the performance of 
the INDO/S method for the case of organic dyes, they 
found that the excitation energies obtained were sig- 
nificantly blue-shifted as compared to the experimental 
results.^ This led the authors to conclude that while the 
INDO/S method performs well for the ultraviolet region 
of the spectrum possibly because it was parameterized 
for small aromatic molecules, its performance is not all 
that good in the visible region of the spectrum.^ There- 
fore, we decided to use the original INDO parameters,^ 
coupled with a high-level correlation scheme (see below) 
so that the influence of parameters is neutralized to a 
certain extent. Indeed, recently El-Shahawy et a/P^also 
used this original INDO approach^ to calculate the exci- 
tation energies of the paracetamol molecule successfully. 
However, to further benchmark our INDO-CI approach, 
as well as to approach the problem from a complemen- 
tary perspective, we have also performed calculations of 
the optical absorption spectra o f the two isomers using 
the ab initio TDDFT approach! 23 1 24 1 This, we believe, 
has helped us tremendously in critically analyzing our 
INDO-CI results. 

The correlated calculations, beyond the INDO-HF, 
are performed using the multi-reference singles-doubles 
configuration-interaction (MRSDCI) approach as imple- 
mented in the computer program MELD. 34 MRSDCI ap- 
proach is a well-established quantum-chemical approach 
in which one considers singly- and doubly-excited config- 
urations from a number of reference configurations, lead- 
ing to a good treatment of electron correlations both for 
the ground and the excited states, in the same calcula- 
tion. Using the ground- and excited-state wave functions 
obtained from the MRSDCI calculations, electric dipole 
matrix elements are computed and subsequently utilized 
to compute the linear absorption spectrum assuming a 
Lorentzian line shape. By analyzing the wave functions 
of the excited states contributing to the peaks of the com- 
puted spectrum obtained from a given calculation, big- 
ger MRSDCI calculations are performed by including a 
larger number of reference states. The choice of the refer- 
ence states to be included in a given calculation is based 
upon the magnitude of their contribution to the CI wave 
function of the excited state (or states) contributing to a 
peak in the spectrum. This procedure is repeated until 
the computed spectrum converges within an acceptable 
tolerance, and all the configurations contributing signifi- 
cantly to various excited states are included in the list of 
the reference states. In the past, we have used such an 
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Figure 1: Structures of: (a) icosahedral (Ih) and, (b) 
quasi-planar (Cs v ) clusters of Bi 2 considered in these 
calculations. 



iterative MRSDCI approach on a number of conjugated 
polymers to perform large-scale correlated calculations of 
their linear and nonlinear optical properties.^ 



III. CALCULATIONS AND RESULTS 



Geometries 



In Figs. [Ta| and [Tb] we depict the structures of the 
icosahedral (Ih) and quasi-planar (Cs v ) isomers of B12, 
respectively. The ground state optimized geometries of 
both the isomers of B i2 were obtained at the INDO-HF 
level. For the icosahedron we obtained the edge length 
to be 1.702 A which is virtually identical to that ob- 
tained by Hayami for the same structure.^ In order to 
further ensure that our geometry was reasonable, using 
the GaussianOSp^ program we also performed geometry 
optimization for both the isomers within the DFT ap- 
proach, employing the B3LYP functional, and a 6-31g(d) 
basis set. This calculation also yielded the edge-length of 
1.70 A for the Ih structure, in perfect agreement with our 
INDO-HF results, and those of Hayami. 11 For the quasi- 
planar structure of Cs v symmetry, we used the geometry 
optimized by Boustani with the bond lengths 1.60 A, 1.64 
A, and 1.65 A. 9 Our own B3LYP/6-31g(d) optimization 
for this system yields corresponding bond lengths to be 
1.63 A, 1.66 A, and 1.68 A, again in good agreement with 
the results of Boustani. 9 However, our results both at the 
INDO-HF and the INDO-CI level predict cage structure 
to be more stable as compared to the planar one, which is 
exactly opposite to the result obtained by BoustaniP In 
a recent experimental study Zhai et al 12 discovered the 
quasi-planar B12 isomer, without ruling out the possibil- 
ity that the larger planar boron clusters will eventually 
fragment into Bi 2 icosahedra. The purpose of our work, 
in any case, is not to examine the relative stability of 
the two isomers, but to investigate their linear optical 
response. 




(d) LUMO 
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Figure 2: (Color online) Molecular orbit als (iso plots) of 
icosahedral B i2 obtained from INDO-HF calculations. 



B. Molecular Orbitals 

Next we present the plots of the molecular orbitals 
(MOs) of both the icosahedral and quasi-planar Bi 2 ob- 
tained from the INDO-HF calculations, in Figs. [2]and[3j 
respectively. Molecular orbitals presented here are close 
to the Fermi level, and corresponding orbitals obtained 
from the first principles DFT/B3LYP calculations are 
presented in Figs. [7] and [5] of the Appendix. Orbital en- 
ergies of some of the INDO-HF orbitals are presented in 
Table [T| Upon comparing the MOs obtained from INDO- 
HF calculations to those of DFT/B3LYP calculations 
presented in these figures, we find perfect agreement as 
to the qualitative nature of the orbitals. Regarding the 
quantitative aspects, the DFT/B3LYP MOs are expect- 
edly more diffuse as compared to the INDO-HF ones, 
because they employ an extended basis set. 

Because of the inversion symmetry of icosahedral B12, 
its MOs are either symmetric (gerade) with respect to 
the inversion operation, or antisymmetric (ungerade) 
with respect to it. From Fig. [2] it is obvious that all 
the MOs from HOMO-2 to LUMO+1 have ungerade 
symmetry while LUMO+2 has gerade symmetry. This 
implies that HOMO^LUMO transition is dipole disal- 
lowed, and, therefore, the closest orbital to which HOMO 
electrons can be optically excited is LUMO+2. Moreover, 
orbitals LUMO and LUMO+1 are degenerate (cf. Table 
|lf because of the symmetry. The charge density distribu- 
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Orbitals 


Orbital Energy (eV) 


Quasi-planar (C?,v) 


Icosahedral (Iu) 


HOMO-5 


-16.227 


-14.749 


HOMO-4 


-15.227 


-14.719 


HOMO-3 


-15.219 


-8.770 


HOMO-2 


-15.219 


-8.731 


HOMO-1 


-10.490 


-8.349 


HOMO 


-10.490 


-3.684 


LUMO 


0.287 


2.475 


LUMO+1 


0.287 


2.475 


LUMO+2 


0.841 


3.476 


LUMO+3 


3.833 


3.593 


LUMO+4 


5.297 


3.633 


LUMO+5 


5.297 


4.028 



Table I: INDO canonical Hartree-Fock orbital energies 
of some of the orbitals close to the Fermi level for both 
the quasi-planar and the cage-like B i2 clusters. 



tion of HOMO (cf Fig. 2c) is concentrated near the two 
pentagons, and vertex atoms on top/bottom of it, with 
negligible charge in between the two pentagons. Thus, 
compared to HOMO, it can be argued that the charge 
distribution of LUMO+2 (cf Fig. [2fJ can be obtained by 
transferring some charge from pentagon and top/bottom 
atoms, to the region between the two pentagons. 

Quasi-planar (C$ v ) isomer of Bi 2 lacks inversion sym- 
metry so that its molecular orbitals cannot be classified 
as gerade or ungerade, a fact which is obvious from the 
MO plots presented in Fig. [3] As a result of this, the 
HOMO^LUMO transition is dipole allowed. However, 
symmetry C$ v of the system manifests itself in form of 
orbital degeneracies, with HOMO being degenerate with 
HOMO-1, and LUMO with LUMO+1 (cf Table 0. The 
charge distribution in HOMO-1 and HOMO is fairly de- 
localized with most regions of the cluster covered, while 
for LUMO and LUMO+1 it is mainly concentrated on 
the edge atoms, with central triangle having negligible 
charge. 



C. Optical Absorption Spectra 

In this section we present the results of our INDO- 
CI and first-principles TDDFT calculations of the linear 
optical absorption spectra of the two isomers. 

Even with a valence electron approximation the num- 
ber of orbitals involved in the INDO-CI calculation is 
rather large (18 occupied and 30 virtual orbitals) and 
can lead to very large CI expansions. We solve this prob- 
lem by freezing occupied orbitals far away from the Fermi 
level. This orbital freezing is carried out in a systematic 
manner, and the convergence of the results with respect 
to the total number of frozen orbitals (Nf reez ) is carefully 
examined in Appendix |B) 



(a) HOMO-1 



(b) HOMO 



(c) LUMO 



(d) LUMO+1 



Figure 3: (Color online) Molecular orbitals (iso plots) of 
quasi-planar B12 obtained from the INDO-HF 
calculations. 



For the cage, we only utilized the inversion symme- 
try (symmetry group C s ), so that the orbitals and the 
many-electron states can be classified into A g (gerade) 
and A u (ungerade) irreducible representations (irreps). 
The ground state belongs to the A g irrep, while the one- 
photon excited states belong to the A u irrep. Our final 
MRSDCI results for the icosahedral structure were ob- 
tained by freezing twelve occupied orbitals (Nf reez = 12) 
both for the A g and A u symmetry manifolds. For the 
quasi-planar isomer, point-group symmetry was not used, 
and thus, the ground, and the excited states, were com- 
puted in the same MRSDCI calculation. The total num- 
ber of CI configurations used in the MRSDCI calculation 
is about one million in case of cage structure, and over 
half a million in case of quasi-planer structure. From 
the sizes of these CI matrices it is obvious that these 
calculations are fairly large-scale, and that the electron- 
correlation effects are properly accounted for. The de- 
tailed analysis of the convergence of our results with re- 
spect to the number of the reference states (N re f) used 
in the MRSDCI calculations, and, therefore, the size of 
the CI matrix, has been provided in Appendix. [B] 

The linear absorption spectrum was computed using 
the sum-over-states approach, with the total number of 
excited states ranging anywhere from fifty to hundred. 
Obtaining that many excited states from large MRSDCI 
calculations was a computationally expensive task, and 
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Table II: Excitation energies, E, and many-particle 
wave functions of the excited states corresponding to 
some of the peaks in the INDO-MRSDCI linear 



absorption spectrum of icosahedral B12 (cf Fig. 4a), 

along with the squares of their dipole coupling 
(/i 2 = J2i \ (f\di\G)\ 2 ) to the ground state. |/) denotes 
the excited state in question, |G), the ground state, and 
di is the z-th Cartesian component of the electric dipole 
operator. In the wave function, the bracketed numbers 
are the CI coefficients of a given electronic 
configuration. Symbols H/L denote HOMO/LUMO 
orbitals. Same information about rest of the peaks can 
be found in table IV of the Appendix [Cj 



Peak 


E (eV) 


/i 2 (a.u.) 


Wave Function 


I 


0.8745 


0.0391 


\H ^L + 5)(0.5818) 








\H -> L + 1]H -> L + 2)(0.4295) 








\H ^L;H ^ L + 3)(0.3355) 








\H ^L;H ^ L + 6)(0.3036) 


II 


1.3642 


0.2215 


\H ^ L + 3)(0.5821) 








\H ^ L + 5)(0.3354) 








\H ^ L + l-H ^ L + 2)(0.3330) 








\H ^ L;H ^ L + 3)(0.3152) 


III 


3.4416 


0.1486 


|if->L + 2) (0.4449) 








\H ^ L + l-H ^ L + 3)(0.3560) 








\H ->L;# ^ L + 2)(0.3496) 








|if->L + 4) (0.3346) 


VI 


6.7030 


0.2619 


|#^L + 10)(0.3260) 








\H ^ L + 1;H ^ L + 8)(0.2449) 








1 -> L + 2)(0.2380) 




6.7366 


0.1417 


|#^L + 9)(0.3778) 








|if -> L; if -> L + 10)(0.2895) 




6.7564 


0.3198 


|if^L + 10)(0.5598) 








\H ^L;H ^ L + 9)(0.3254) 


VIII 


7.7822 


0.3373 


|if-3^L + 6)(0.3285) 








|if-3->L + l;iJ->L + 4)(0.2601) 








\H-3^L;H ^ L + 3)(0.2586) 




7.9592 


0.1889 


|i?->L + l;#-3->L + 6)(0.2851) 








|if-2^L + 2)(0.2586) 








|if - 1 -> L;# -> L + 1; 
if->L + 4)(0.2514) 




8.0153 


0.1343 


\H-3^L;H ^ L + 6)(0.2837) 








|if-3^L + 5)(0.2096) 



took several days for some calculations on our computer 
cluster, running dual-CPU quad-core processors. Our 
INDO-MRSDCI linear absorption spectra of the cage and 
planar isomers are presented in Figs. [4] and |5j respec- 
tively. With the purpose of benchmarking our INDO-CI 
approach, and also to study these systems from a com- 
plementary perspective, we performed ah initio TDDFT 
calculations of the absorption spectrum of both the iso- 
mers. For the purpose we used the same geometry as 



Table III: This table contains information pertinent to 
some of the peaks of INDO-MRSDCI optical absorption 
spectrum of the quasi-planar Bi 2 as shown in Fig. [5a] 
The symbols have the same meaning as in the caption 
of table [TTJ Same information about rest of the peaks 
can be found in table [V] of the Appendix [Cj 



Peak 


E(eV) 


/i 2 (a.u.) 


Wave function 


I 


4.6844 


0.0960 


|if-l->L + l)(0.7525) 








|if->L)(0.5714) 








\H - 1 -> L; H -> L + 1)(0.1271) 


II 


8.4488 


0.3599 


\H -> L + 3)(0.8665) 








\H-3 L)(0.1808) 








\H -> L + 4)(0.1673) 








\H- 1 L + 5)(0.1361) 








|if -4->L + l)(0.1278) 








\H L + l-H L + 2)(0.3330) 


III 


8.9674 


0.1690 


1 -> L + 2)(0.4386) 








1 -h. L + 5)(0.4044) 








|if -> L + 4)(0.3889) 








|if -3^L + 2)(0.3457) 








|if-2->L + l)(0.2625) 








|#-4^L + 1)(0.2528) 








\H ^ L + 6)(0.2008) 


IV 


9.4482 


0.0505 


|#-3^L + 1)(0.6229) 








|if-2->L)(0.4281) 








|if - 1 -h. L; if -h. L + 1)(0.2539) 








|if-l->L + 4) (0.2457) 








|if -2^L + 2)(0.2268) 



those for the INDO calculations, and used the GAUS- 
SIAN 03 package, 39 coupled with the 6-31g(d) basis set, 
and B3LYP gradient hybrid correlation functional. These 
TDDFT absorption spectra are also presented in Figs. [4] 
and[U 

Upon comparing the INDO-MRSDCI results with the 
ah initio TDDFT ones, we conclude that the spectra 
computed by the two approaches are in very good qualita- 
tive agreement with each other, with both sets of spectra 
exhibiting weak absorption at lower energies, and very 
intense absorptions at high energies. Next, we make a 
quantitative comparison between the spectra computed 
by the INDO-MRSDCI and TDDFT approaches, for 
both the isomers. 

For the icosahedral cluster (cf Fig. [I]) the two spectra 
exhibit the onset of optical absorption close to 0.9 eV in 
form of low intensity peaks, (b) next set of low intensity 
peaks starts a little below 4 eV in both the spectra; in the 
INDO spectrum these peaks continue beyond 4 eV, while 
in the TDDFT spectrum they are all below 4 eV, and (c) 
high energy feature of both the spectra are dominated 
by very high-intensity peaks which occur above 6 eV in 
the INDO spectrum and between 4 eV and 5 eV in the 
TDDFT spectrum. Thus, we have very good agreement 
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(a) INDO-MRSDCI spectrum 
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between the spectra for the first set of peaks around 0.9 
eV and the second set of peaks at energies close to 4 
eV. For high-energy features we have a disagreement in 
that the INDO theory predicts these features at energies 
higher than 6 eV while the ah initio TDDFT predicts 
them between 4 and 5 eV. 

Upon comparing the INDO-MRSDCI and TDDFT 
spectra for the quasi-planar isomer (cf Fig. [5]), a similar 
picture emerges. In the INDO-MRSDCI spectrum the 
first, relatively weak, peak occurs at around 4.7 eV fol- 
lowed by a set of strong peaks beyond 8 eV, with the inter- 
mediate region (4.7 eV — 8.4 eV) exhibiting almost negli- 
gible absorption. The TDDFT spectrum is slightly differ- 
ent with the first weak peak close to 4.04 eV, followed by a 
slightly stronger peak at 5.27 eV. Really intense peaks in 
the TDDFT spectrum occur in the energy range between 
5.9 eV and 7 eV. Therefore, the quantitative comparison 
between the INDO-MRSDCI and TDDFT spectra for the 
quasi-planar isomer is quite similar to that of the icosahe- 
dral isomer. Agreement is reasonable for the lower energy 
peaks, but for the higher peaks, TDDFT spectrum is sig- 
nificantly redshifted as compared to the INDO-MRSDCI 
spectrum. Qualitatively, however, both sets of spectra 
are in very good agreement with each other. 

Next we present a detailed comparative analysis of our 
INDO-MRSDCI results for the cage and the quasi-planar 
structures. To facilitate this comparison visually, we also 
present a combined plot of the two spectra in Fig. [6] 
A cursory look at Fig. [6] reveals that lower energy re- 
gions of the absorption spectra consist of relatively lower 
intensity peaks while the higher energy spectra in both 
the isomers exhibit intense absorption. A comparison of 
the the spectra of the two isomers reveals that: (a) the 
absorption spectrum of the icosahedral isomer is signifi- 
cantly red-shifted as compared to the planar one, and (b) 
the peak intensities in both the spectra are of the same 
order of magnitude. The main distinguishing feature of 
the absorption spectra of the two isomers, as depicted by 
our calculations, is that in the icosahedral B12 the optical 
absorption (cf. Fig. [3 begins at rather low energies with 
the first peak (peak I) located at 0.86 eV, with two more 
peaks (peaks II and III) located below 4 eV, while in the 
quasi-planar isomer (Fig. [5| the first absorption absorp- 
tion feature (peak I) is located significantly higher at 4.73 
eV. This INDO-MRSDCI result is also confirmed in the 
TDDFT calculations, which also predict the absorption 
to commence at much lower energies in the icosahedral 
isomer, as compared to the quasi-planar one (cf. Figs. 
4b and 5b). This prediction of our calculations can be 
tested in future experiments on boron clusters, and if 
confirmed, can be used for distinguishing the icosahedral 
clusters from the quasi-planar ones in optical absorption 
experiments. 

As far as the higher energy features are concerned, in 
the icosahedral cluster (Fig. [4|, a series of high intensity 
features (peaks V, VI, VII, and VIII) start around 6.7 
eV and continue past 8 eV. In the planar isomer (Fig. [5| 
on the other hand the higher energy absorptions (peaks 



(b) TDDFT spectrum 



Figure 4: Linear optical absorption spectrum of 
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II through VII) start beyond 8 eV and continue to much 
higher energies beyond 12 eV. 

The many-particle wave functions of the excited states 
corresponding to some of the peaks in the spectra, along 
with the squares of their transition dipoles, are presented 
in Tables [TT] and IIIII The same information about the 
rest of the peaks can be examined in tables [IV] and [V] 
of Appendix [Cj One distinguishing feature of the optical 
absorption in cage, compared to that in the quasi-planar 
isomer is that HOMO (H) to LUMO (L) transition is 
dipole forbidden in the cage because, as discussed above, 
H and L orbitals have the same inversion symmetry of u, 
and among the unoccupied orbitals, L + 2 is the lowest- 
lying with the opposite symmetry of g, thereby making 
H — > L + 2 as the lowest allowed single-particle optical 
transition for the cage. For the quasi-planar B12 on the 
other hand, the single-particle transition H — > L is opti- 
cally allowed because the system does not have inversion 
symmetry. 

Before further discussions of our results it is useful to 
remember that in the crystalline form boron is a semi- 
conductor, with an indirect band gap close to 1.53 eV 
for /3-rhombohedral boron, 36 whose crystal structure is 
also based on Bi 2 icosahedron. Semiconductors exhibit 
single-particle inter-band optical absorption at lower en- 
ergies, and collective plasmon absorption at high ener- 
gies. Therefore, it is of interest to understand the op- 
tical absorption spectra of B i2 clusters based on these 
pictures and to ascertain whether absorption exhibits 
single-particle behavior or the collective one. Indeed for 
the case of metallic clusters such as those of alkali met- 
als the issue, whether the optical excitations are similar 
to those in a molecule, or they are related to the plas- 
mons in the bulk, has been extensively examined!^ In 
metal cluster physics plasmon excitations are identified 
by appearance of continuous absorption pattern. More- 
over, Koutecky and coworkers 38 have devised a scheme 
according to which if the many particle wave function of 
a given excited state is dominated by one singly-excited 
configuration it is classified as a normal inter-band ab- 
sorption. 

On the other hand if the excited state wave func- 
tion is a linear combination of several configurations 
with similar weights, it is called a collective (plasmon 
like) excitation. 38 Close examination of the many-particle 
wave functions of both the isomers presented in Tables 
PH |III| |IV[ and [V] reveals that none of the excited states 
participating in the optical absorption of B12 clusters 
are dominated by single configurations. All these ex- 



cited states have prominent multi-reference character, 
and are linear combinations of several singly- and doubly- 
excited configurations with significant coefficients. For 
example, peaks I of both the isomers are a mixture of 
singly and doubly excited configurations, as are all the 
other peaks. This, as per the criterion outlined in the 
work of Koutecky and coworkers, 38 suggests plasmon 
like behavior. Moreover, the absorption spectrum ap- 
pears to be fairly continuous starting with low-energy 
and low-intensity absorptions, followed by high-energy 
high-intensity absorptions before tapering off. 



IV. CONCLUSIONS AND FUTURE 
DIRECTIONS 



In conclusion, we have presented a theoretical study 
of linear optical absorption in two B12 isomers, namely 
icosahedral and quasi-planar clusters. The wave- 
function-based calculations were performed using the 
semi-empirical INDO model pi and electron correlation 
effects were taken into account by means of large-scale 
MRSDCI computations. To obtain an alternative per- 
spective on the optical absorption, TDDFT method 
based ab initio calculations of the spectrum were also 
performed. A comparison of the INDO-MRSDCI and 
TDDFT calculations revealed that the two methods lead 
to spectra which are qualitatively very similar. On the 
quantitative front, the high-energy features of spectra 
were found to be red-shifted in the TDDFT approach 
as compared to the INDO-MRSDCI approach. Which 
of these results is closer to reality can only be decided 
by the experiments, which, we hope will be performed in 
the future. Another aspect of these spectra is that the 
optical absorption in the icosahedral clusters begins at 
much lower energies as compared to the planar one, a fact 
which can be used in the optical detection of these clus- 
ters. The high-intensity absorption in the cluster takes 
place at higher energies, which our calculations suggest, 
are plasmonic in nature. Results of these calculations, 
which to the best of our knowledge have not been per- 
formed earlier, can be tested in the future experiments. It 
will also be interesting to investigate the nature of triplet 
excited states in boron based clusters along with their 
nonlinear optical properties. Calculations along those 
directions are in progress in our group, and results will 
be presented in future publications. 
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Appendix A: MOLEULAR ORBITALS OF B 12 
ISOMERS OBTAINED FROM AB INITIO 
CALCULATIONS 

In this section we present some of the molecular or- 
bitals of icosahedral and quasi-planar isomers of B12, ob- 
tained from the first-principles DFT calculations employ- 
ing B3LYP functional, and the 6-31g(d) gaussian basis 
set. The calculations were performed using the Gaus- 
sian03 program,^ and the orbitals presented are the ones 
close to the Fermi level. In Fig. [7] orbitals of icosahedral 
isomer are plotted, while in Fig. [8] those of the quasi- 
planar isomer are presented. These orbitals agree quite 
well with the corresponding orbitals obtained from the 
INDO-HF calculations, presented in Figs. [2] and [3] of the 
main text. 



Appendix B: CONVERGENCE ISSUES 

In this section we discuss the convergence of our MRS- 
DCI calculations with respect to: (a) the number of 
frozen orbitals Nf reezj and (b) the number of reference 
configurations N re f, with respect to which the singly- 
and doubly-excited configurations spanning the CI space 
are generated. For a given value of Nf reez , successively 
larger MRSDCI calculations (i.e. with the increasing val- 
ues of N re f) were performed, until no significant changes 
in the results were observed both with respect to, decreas- 
ing Nf reez , and increasing N re f. In Fig. [9] we present the 
results of the best (largest N re f) MRSDCI calculations 
on the icosahedral B12 with Nf reez = 16, 14, and 12, and 
N re f = 83, 92, and 84, respectively. In these calcula- 
tions the sizes of the MRSDCI matrices (N total) ranged 
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from 20935 for N freez = 16 to 940945 for N freez = 12. 
From the figure it is obvious that there are no signifi- 
cant qualitative or quantitative changes in the absorption 
spectrum as Nf rP ^ z is decreased from fourteen to twelve. 



Similarly, Fig. [10] demonstrates the convergence of the 
absorption spectrum with respect to increasing N re f (and 
hence N tot ), with values N re f = 1, 60, 79, and 84. Thus, 
we can conclude that our results for the cage B12 are well 
converged both with respect to Nf reez as well as N re f. 

Similar checks of convergence were also performed for 
the quasi-planar isomer of B12, figures corresponding to 
which are not presented here for the reasons of brevity. 



Appendix C: EXCITED STATE WAVE 
FUNCTIONS, ENERGIES, AND TRANSITION 
MOMENTS 



In the following two tables we present the excitation 
energies, the many-particle wave functions, and the tran- 
sition dipole moments with respect to the ground state, 
of the excited states corresponding to those peaks in the 
INDO-MRSDCI linear absorption spectra of the two iso- 
mers of B12, which were not included in the tables [TT] and 

urn of sec. nni 
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(a) INDO-MRSDCI spectrum 




(b) TDDFT spectrum 

Figure 5: Linear optical absorption spectrum of 

quasi-planner Bi 2 computed using: (a) the 
INDO-MRSDCI method, with N freez = 12 and 
N ref = 31, and (b) TDDFT method using 6-31g(d) 
basis set and B3LYP functional. A line width of 0.1 eV 
was used to compute the spectra in both the cases. 
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(eV) 



Figure 6: (color online) Comparison of the linear optical 
spectra of icosahedral (black) and quasi-planar (red) 
structures of I^cluster. 



12 




Figure 7: (color online) Some molecular orbitals (iso 
plots) of icosahedral Bi 2 , obtained from the first 
principles DFT/B3LYP calculations. 
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(a) HOMO-1 (b) HOMO (c) LUMO 



(d) LUMO+1 



Figure 8: (color online) Some molecular orbit als (iso 
plots) of quasi-planar Bi 2 , obtained from the first 
principles DFT/B3LYP calculations. 
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Figure 9: (color online) Convergence of the linear 
absorption spectrum of icosahedral B12 computed using 
the INDO-MRSDCI method for decreasing number of 
frozen occupied orbitals. A line width of 0.1 eV was 
used to compute the spectra. 



15 




Figure 10: (color online) Convergence of the linear 
absorption spectrum of icosahedral B12 with increasing 

number of reference configurations (N re f) in the 
MRSDCI wave function with Nf reez = 12. A line width 
of 0.1 eV was used to compute the spectra. 
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Table IV: Excitation energies, E, and many-particle 
wave functions of the excited states corresponding to 
some of the peaks in the INDO-MRSDCI linear 



absorption spectrum of icosahedral B12 (cf. Fig. 4a), 

along with the squares of their dipole coupling 
(/i 2 = J2i \ (f\di\G)\ 2 ) to the ground state. |/) denotes 
the excited state in question, |G), the ground state, and 
di is the z-th Cartesian component of the electric dipole 
operator. In the wave function, the bracketed numbers 
are the CI coefficients of a given electronic 
configuration. Symbols H/L denote HOMO/LUMO 
orbit als. 



Peak 


E (eV) 


/i 2 (a.u.) 


Wave Function 


IV 


4.1563 


0.0919 


\H — > L 4- 3U0 4902) 








\H — > L;H^ L + 6)(0.4622) 








\H ^L;H ^ L + 5)(0.4070) 


V 


5.0252 


0.0171 


\H-3^L;H ^ L + 2)(0.3519) 








\H - 2 L]H -> L + 2)(0.3322) 


VII 


7.0917 


0.3664 


\H ^ L + 10)(0.5795) 








\H ^L;H ^ L + 9)(0.4872) 








|if->L + l;if->L + 9)(0.3962) 


IX 


8.2770 


0.1232 


\H-2^L;H ^ L + 5)(0.3093) 








\H-2^L;H -h. L + 3)(0.2999) 




8.3141 


0.0905 


|if->L + l;if-l->L + 2)(0.3300) 








|if-l->L + l;if->L + 4)(0.2080) 








1 -> L + 3)(0.2043) 




8.3178 


0.0974 


\H -2^ L + l-H ^ L + 2)(0.2716) 








\H -> L + 4; if - 1 -> L)(0.2525) 








|if-l^L + l;if^L + 65)(0.2346) 








\H -2-> E,H -> L + 4)(0.2121) 
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Table V: This table contains information pertinent to 

some of the peaks of the INDO-MRSDCI optical 
absorption spectrum of the quasi-planar Bi 2 , as shown 
in Fig. [5a| The symbols have the same meaning as in 



the caption of table IV 



Peak 


E(eV) 


/x 2 (a.u.) 


Wave function 


V 


10.7578 


0.0279 


\H ^ L + 8)(0.6230) 








|if-5^L + 2)(0.4956) 








\H L + 1; if -> L + 1)(0.2400) 








\H - 5 L)(0.1807) 


VI 


11.2054 


0.2187 


\H - 1 -> L + 9)(0.4144) 








|if -5 -h. L)(0.3539) 








|if -> L + 10)(0.3255) 








|if- 5 ^L + 2) (0.3049) 








lif ^ L + 8)(0.2973) 

1 1 / V / 








|if -> L;iJ -> L)(0.1953) 


VII 


12.3562 


0.1598 


Iff -> L + 10)(0.6469) 








|if-3^L + 3) (0.3408) 

1 1 / V / 








|if -> L + 12)(0.2985) 








|if -> L + 13)(0.2299) 








|#^L + 8)(0.2095) 








|iJ-l->L + ll)(0.2063) 








\H -4^ L + 5)(0.1792) 


VIII 


13.1344 


0.0300 


|if-3^L + 3)(0.6665) 








|if-4-> L + 5)(0.3184) 








|if -> L + 13)(0.2582) 








|tf-2^L + 2)(0.2586) 








\H - 1 -h. L; H -> L + 3)(0.2186) 








|if-4-> L + 7)(0.1921) 








|if - 4 L; if -> L + 2)(0.1574) 



